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FINAL REPORT NASA Grant NAG-1-757 


The final phase of modeling solar pumped lasers considers the photodissociation of 
the perfloralkl molecules n — C 3 F 7 I ,t — C 4 FqI and i — C 3 F 7 I. Computer modeling has 
been compared to laboratory data and good agreement between theory and experiment 
has been achieved. The following is a summary of the modeling of solar pumped lasers. 

A perfloride gas enters a tube at a point z = 0 and travels a distance L at a velocity 
w and then exits at the point z — L. During the flow the gas is solar pumped over the 
initial distance 0 < z < z 0 , where z 0 is less than or equal to L. The perfloride gas interacts 
chemically with light and a chain of chemical reactions occur. 


CHEMICAL KINETICS 

Following reference 1,2 and 3, let RI denote one of the perfluoralkl molecules, where 
R stands for either n — C 3 F 7 or t — C 4 Fq. The perfloralkl molecule flows through a circular 
tube where the gas is radially illuminated and pumped by the ultraviolet light in the 270 
nm range and is characterized by the chemical reaction 


RI + hu 


270 nm 


R + T , 


where I* is the 2 Pi / 2 excited state of iodine, which emits light at 1.315^ m. In addition, 
to the photodissociation reaction there are other chemical reactions which occur. These 
reactions are summarized as follows. 




1 



Photochemical Reactions for Iodine Laser Simulation 
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Stimulated emission 


Regeneration reactions 


(rates have units [cm 3 / sec]) 
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Quenching reactions (rates have units [cm 3 /scc]) 

7* +72/ I + RI 

r + / 2 i + h 

I* + R -&* I + R 
I*+R 2 / + / 2 2 

r + / / + / 


( 1 ) 


( 2 ) 


(3) 
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Three body collision reactions (rates have units [cm 6 /sec]) 


r + i + Ri 

i 2 + ri 

I + I + RI -2* 

i 2 + ri 

r + 1 + 12 

i 2 + 1 2 

i + i + h 

i 2 + 1 2 

i + i + r 2 

I 2 + r 2 

i + i* + r 2 

♦ i 2 + r 2 


Wall reactions (rates have units [cm 3 jsec]) 

I + I + Wall I 2 + Wall 
I+R 2 +Wall R + RI + Wall 

Pyrolysis reactions (rates have units [cm 3 /$ec]) 


RI+ -^R + I 

R 2 ^R + R 

In the above reactions note that we have denoted by 
g,, i = 1, ...9 those reaction rates which are quenching; 
ki,i = those reaction rates which are regenerative; 

ki, i = 9, 10 those reaction rates which are due to pyrolysis; 
v,-,* = 1,2 the reaction rates of wall reactions; 

Ci, i = 1,6 those reaction rates due to three body collisions. 

The dimension of the reaction rates are s -1 for single particle decomposition; 
for binary collisions and cm 6 s -1 for three body collisions. 


RATE COEFFICIENTS 


The reaction rate coefficients are assigned the following values for n — C 3 F 7 I. 
Regenerative and pyrolysis 


ki = 1.0 (10) -14 

fc 6 = 0.0 

k 2 = 2.3 (lO) -11 

k 7 = 3.0 (10) _19 

k 3 = 6.5 (10) _13 

= 1.6 (10) ~ 23 

k 4 - 3.0 (10) ~ 16 

kg = 1.0 (10) 15 exp(— 2.48((10) 4 ) /T) 

k 5 = 5.0(10)- n 

k 10 = 1.0 (10) 17 exp(— 4.73((10) 4 ) /T) 


Three body 

ci = 1.6 (10)" 33 

c 2 = 5.7 (10) -33 exp(1360/r) 
c 3 = 0.0 

c 4 = ( io)- 29 - 437 - 6 - 844io » lo ( T / 3OO )+ 2 - 163 ( L<> 0 lo ( r / 300 )) 3 
c 6 = 8.0 (10) " 33 exp(1310/T) 

Cq — 0.0 

Quenching 

gl = 1.7 (10)" 17 

q2 = 2.9 (10)" 11 exp(-4.4(10)" 3 (T - 300)) 

<73 = 3.7 (10) -18 
q4 = 4.7 (10) _ 16 
q5 = 1.6 (10)~ 14 
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Pumping and Lasing 


rl = 0.0 


v2 = 0.0 

^ = 1.2(lO)- 2 (n-C 3 F 7 I) 

£2 = 0.12 (n - C 3 F 7 I ) 

6 = 0.0 


Note that many of the coefficients are temperature dependent, where T is the absolute 
temperature. In addition, certain reactions, although postulated, are not readily known. 
For example, the effect of wall reactions. These postulated reactions, though used in some 
of the preliminary studies, have been set to zero by setting the reaction rate coefficients 
to zero. The £1 and £2 coefficients are multiplied by the pump power density in solar 
constants and £3 = 0 coefficient is associated with one of the postulated reactions which 
was removed. The above reaction kinetics are in basic agreement with those given in 
references 1,2, 3, 4 and 5. 


DIFFERENTIAL EQUATIONS FOR REACTION KINETICS 

Using the law of mass action, the above set of chemical kinetics are replaced by a 
coupled system of ordinary differential equations. We denote the concentration of the 
various chemical reactants by 

xi = [RJ], x 2 = [R], x 3 = [R 2 ], x 4 = [J 2 ], 

X 5 = [/*], X 6 = [/], X 7 = p+, x 8 = P- 
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where p+ and p- are photon densities to be discussed shortly. The chemical kinetics can 
now be represented in the form 


dxi 

dt 


^ + (V • V) i,- + x t V . V = /,•(*), * = 


. . — * 

where, £ = (zi, x 2 , X3, x 4 , ig, x 6 ) and V = u e z is the flow velocity in the z direction. Here 
(1? • V)x, is a convective term and x,- V -V is a dilatation term. For steady state conditions 
we write the system in the form 


dxi . . , doj . , , 

w — = fi(x) - Xi— = F,(x), * = 1, ...,6. 


where = 1,...,6 represents the right hand side of the above system of equations and 
are given by 
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fi(x) =fcx[J?][J*] + * 2 [*M - k 7 [r\[RI] - k 4 \R][RI } 

- ke[R][RI] - hlRI} + v 2 [R 2 ][I) + k s \R){I 2 ] - k 9 [RI ] 
h(x) =6 - fci[i2][r] - k 2 [R][I] - 2k z [R} 2 - k 4 [RI][R] 

- Jfc 6 [727][72] - ksMh] + v 2 [R 2 ][I] + k 7 [RI]\r] + fc 8 [/][i 2 /] + k g [RI] 

+ 2fcio [-^ 2 ] — [72]/l" 2 

f z {x) =k z [R ] 2 + ke[RI][R] + k 4 [RI}[R] - v 2 [R 2 ][I] - M[72 2 ] - [R 2 }/t z 
A(x) = - hlh] + CtlRIlimi] + c 2 [RI][I ] 2 + c 3 [I 2 ][I*][I] + c 4 [7 2 ][7] 2 
+ k 7 [RI][I*} — A: 5 [iE][/ 2 ] + vi[I} 2 
+ c 5 [7] 2 [72 2 ] + k 8 [RI][I] + c 6 [7][7*][72] - [I 2 \/t 4 

f*(x) =Qytl[Rl} + 0.5ie 2 [/ 2 ] - *1 mn - ? 2 [/ 2 ][/*] - cc([p + ] + [/»-])([/*] - \{I)) 

+ ke\R}[RI} - g 3 [72][7*] - ? 4 [72 2 ][7*] - <7s[7*][7] - k 7 [RI][I*} 

- c 6 [iE 2 ][r][7] - cxMrn/] - c z \i 2 ][i*][i) - qi min - [r]/r 5 
fe(x) =1.496[/ 2 ] + ?i[727][7*] + ? 2 [/ 2 ][/*] - 2c 6 [7] 2 [72 2 ] - k a [I][RI] 

cc([p + ] + [„_])([/*] - \[I}) - ci[727][7*][7]-2c 2 [727][7] 2 - c 3 [J 2 ][/*][J] 

- 2c 4 [7 2 ][7] 2 - fc 2 [i?][/] + k 4 \RI][R] + q z [I* ][72] + g 4 [7*][72 2 ] 

+ g 5 [r][/] + k 5 [R]\I 2 ] - v 2 [R 2 }[I] - 2 t>i[/] 2 - c 6 [72 2 ][7*][7] + k 9 \RI } - [/]/r 6 
/ 7 (x)=ca[p + ]([r]-i[/]) 

(7) 

where Q y is the quantum yield and the terms involving r,, » = 2,..., 6 represent diffusion 
effects. Without the diffusion terms, the above system of nonlinear differential equations 
conserves the masses of the species involved in the reactions and for steady state operation 
at any point z we have the immediate integrals 

[727] + [72] + 2[72 2 ] =Constant 

( 8 ) 

[727] + 2[7 2 ] + [7*] + [7] =Constant. 
provided the diffusion terms are removed. 

The above coupled system of nonlinear ordinary differential equations are subject to 
the initial conditions that at z = 0 we have 

xi = B, X{ — 0.0 for » = 2,...,6, 
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where B = 9.66 (10) 18 P/T is the initial concentration of [i2/] which is dependent upon 
the inlet pressure P (torr) and absolute temperature T. 

The above system of differential equations is scaled in both the dependent and inde- 
pendent variables. The steady state independent variable is z which varies between 0 and 
L. This distance is scaled to the interval 0 < s < 1 by making the substitution z = sL. 
The dependent variables are scaled by the substitutions 

Xi = Byi, t' = l,... 7, 

and the initial conditions upon y,-,t = 1,...,6 become yi(0) = 1.0 and y,(0) = 0.0 for 
* = 1 , ..., 6 . 

In addition to these chemical variables let p denote the photon density in the laser 
tube which has a length L. Assume p = p_ + p + , where p + is the number of photons per 
cubic centimeter moving to the right and p_ is the number of photons moving to the left. 
From reference 1, the photon density satisfies the differential equation 

—, + ([r]- lyi) 

If ^ ( |r| - 5 |z| ) 

where c is the speed of light, a is a capture cross section, [/*] is the concentration of excited 
state of iodine, and [/] is the concentration of iodine. 

Photon Density Equation 


Let the light flux density of the lasing photon be denoted by p = p + + p- where 
p+ = p+(z,t) denotes the photon density propagating in the positive z direction and 
P- = p-(z,t) denotes the photon density propagating in the negative z direction. The 
differential equations for these photon densities are 


1 dp+ dp+ 
c dt dz 
1 dp- dp- 

c dt dz 


= ffM[r] - i|/|) 

= CT# >_([/-] - ![/]) 


where c is the speed of light in the optical medium, and a is a capture cross section. In the 
above equations a[I*\p + is the amplification factor resulting from population of the upper 
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laser level of the active medium and — |<r[/]p+ is the decrease in photon density due to 
population of the lower lasing level. 

Note that for steady state conditions the quantities p+ and p- satisfy the differential 
relation 


dp- , dp+ _ n 

l ’*~dr +p -~dr-° 

and consequently by direct integration we find that 


(9) 


p+{z)p-[z) = Ko = Constant. (10) 

Denote by Ri and i? 2 the reflectivities of the mirrors then we will have the boundary 
conditions 

p + (0) = .Rip_(0) and P-{L) = R 2 p+{L), 
as illustrated in the figure 1. 



Figure 1. Geometry for photon density boundary conditions 


Therefore, we require that at z = 0 the boundary conditions 

/j + (0)/j_( 0) = i?ipl(0) = -^4(0) = K 0 . 

K\ 

and at z = L we have the boundary condition 

P+(L)p.{L) = = *2f>5-(i) = K 0 . 

k 2 


( 11 ) 


( 12 ) 
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These conditions require that 


P+(0) = y/K^R[ and p + {L) = (13) 

where Ri and R 2 are the mirror reflectivities. 

For steady state operation the initial photon density at z = 0 is unknown. Therefore, 
a shooting method is used to determine this boundary condition. That is, we first scale 
the above equations and introduce the variables for the density p + and 3/7 for its scaled 
value. We then assume a value for the photon density, X7 (or 3/7), at z = 0 and then solve 
the scaled system of differential equations over the interval 0 to 1 and test the value of the 
photon density at the end of the laser tube by using the following arguments. 

For fixed values of the reflectivities we guess at an initial value for K 0 and then 
integrate the system of differential equations to obtain a calculated value for p+{L) which 
is then compared with the theoretical value from equation (13). If the calculated and 
theoretical values are different we iterate on Kq until the final value of P+[L) agrees with 
the theoretical value. This describes the shooting method referred to earlier. Once Ko is 
determined we can use the value of the transmission coefficient given by t m = 1 — R 2 and 
find the output power from the relation 

P = e^tmcp+iL) ( Wjcm 2 ). (14) 

where is the radiation energy quantum from I* given by = 1.5 (10)“ 19 W s. 


Two Pass Amplifier 


The above ideas can also be extended to model a two-pass amplifier. The geometry 
for the two-pass amplifier is illustrated in the figure 2. 
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Figure 2. Geometry for Two-Pass Amplifier. 


In this figure P, n is the power density into the amplifier [Watts /cm 2 }, W is the 
transmission coefficient, e u is the radiation energy from I* [ Watt — sec], c is the speed of 
light, P out is the output power density, Rl is the reflection coefficient of the mirror, p + (z) 
is the photon density for wave motion in the positive direction [cm -3 ], p~[z) is the photon 
density for wave motion in the negative direction [cm -3 ]. 

The input power density P* n is assumed to be known. Also the photon density satisfies 
for all values of z the equation (10) where Kq is a constant. At z = 0 we have 

P+W/M®) = K 0 (15) 

and at z — L we have 

p + (L)p.{L) = K 0 . (16) 

The value for p+(0) is determined from the input power density P, n by the relation 

MO) = ^ (17) 
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(18) 


and consequently from the relation (15) we have 


Mo) 


K o 

/>+(»)' 


In the above equations W is the transmission coefficient of the Brewester windows and Kq 
is an assumed initial value. At z = L we have R% = W 2 Rl as the reflection coefficient 
from the Brewster window. We then are able to calculate the return photon density from 
the Brewster window at z = L. We find this density has the value 


p_(L)=R 2 p + (L). (19) 

Combining the equations (15), (16) and (19) we determine that 

RtfX(L)~Ko or 0+W = = \] < 2 °) 

Using the value given by equation (20) as a theoretical value for p + (L) we can compare this 
value with the calculate value for p+(L) obtained from an integration of the differential 
equations defining the chemical reactions occurring in the amplifier. We adjust the value 
for p(0) (by adjusting Ko) until 

\p+ {P) theoretical P+i.^) calculated] ^ E (2l) 

where E is an error criteria. The above procedure is again a shooting method for the 
determination of p- (0) . When the error criteria is satisfied, the output power is obtainable 
from the relation 

Pout = p ( 0)We„C ( 22 ) 

and the gain is calculated from the relation 

gain = 10 Log 10 [decibels]. (23) 

Again the length of the amplifier is normalized to unity and the system of differential 
equations are solved using a Runge-Kutta integration scheme. Previous progress reports 
have the computer codes for both the iodine laser and the amplifier models. 

In both the iodine and amplifier models the compressibility of the perflouride gas 
requires additional considerations. If the perflouroalkyl iodide gas is incompressible the 
above equations hold. For a compressible gas we must add to the above equations an 
equation of state, a continuity equation, a momentum equation and an energy equation. 
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COMPRESSIBILITY EFFECTS 


Equation of State 

We assume an equation of state for an ideal gas and write 

P = rjRT (15) 

where P is the gas pressure [Pa], »7 is the gas density [Kg/m 3 ], T is the absolute tem- 
perature \K\ , and R = R*/M [J/KgK] with M the mass of the gas and R* is the gas 
constant, R* = 8.317 [J/mole K\. For n — C 3 P 7 / we use the value M = 296 Kg /mole and 
for t — C 4 FqI we use the value M = 346 Kg /mole to convert R* to R. 

Continuity Equation 

The continuity equation expressing the conservation of mass flow can be represented 

^ + div (rjV) = 0 (16) 

where rj is the fluid density [Kg/m 3 ], and V is the fluid velocity [m/s]. For steady state 
conditions and a flow in the axial direction we let V — u)k and reduce the continuity 
equation to the form 

A(„ w ) = o. 

An integration of this equation gives 

r]u = Ci = Constant, (17) 

where C\ is a constant of integration. 
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Momentum Equation 

The momentum equation for a control volume having a mass rjdr, where dr is an 
element of volume, is given by 

M = jjjvndr. 

Using Newton’s second law we write 


F = 


DM 

Dt 


MI 


Vrj dr 


(18) 


where, Dj Dt is the material derivative. Consequently, we find 


DM 

Dt 


= If + Vl dr 


(19) 


where da is a vector element of surface area. The surface integral term in equation (18) 
represents the efflux of momentum through the control volume and the volume integral 
term in equation (18) represents the change in momentum inside the control volume. Using 
the Gauss divergence theorem, we can change the surface integral to a volume integral and 
write 

JJ V{r)Vdo) = JJJ [v-ijfujdr. (20) 

Then the momentum equation becomes 


F = 


DM 

Dt 


= HI AbrfavV) ir + ^ /TYfrw i) 


( 21 ) 


—* —f A A ^ A 

where VV is the dyadic ur kk of V = a> k> and 


V-r,VV = ~{r,ij 2 )k. 


( 22 ) 
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We also make the substitutions F = /// f dr where / is the average force per unit 
volume as this allows the momentum equation to be written in the form 


= III fiT = £ 5T = lil [s ( ’ ,u,) ' k + lz { ^ ] * 


dr 


(23) 


Neglecting the viscosity and using f = —VP as the average force per unit volume (which 
is due to the fluid pressure P), we have 


dP d . . d , 2 \ 


(24) 


From the equation (17) we use rjoj = C\ and examine the steady state form of the equation 
(22) to obtain 

a p a,., 

(25) 


dP _ _ du> 


Integrating the equation (25) produces the result 


P + C 1 u = C 2 , 


(26) 


where C 2 is a constant of integration. 

Energy Equation 

In terms of the specific enthalpy h per unit mass, the energy equation for the fluid 
flow is 


Dh DP 

V-ftt = 


+ /cV 2 r + g 


(27) 


Dt Dt 

where P is the pressure, T is the absolute temperature, k is the thermal conductivity and 
q = q(z) is the radiation heat flux. In one dimension the equation (25) reduces to 


dh 


dh dP 


dP d 2 T 


’ , ar + ’" J aI = '5F +u a7 + ' c ^ 


+ 9- 


(28) 


The specific enthalpy h can be expressed in terms of the specific heat at constant 
pressure and 

(29) 

c; — C* = R* (30) 


h = C p T. 


Using the relation 


15 


the specific heat at constant pressure C* can be expressed in terms of the specific heat at 
constant volume C * and the universal gas constant R * , where R* = 8.317 [ J/moleK ]. We 
use the following numerical data to obtain a least squares approximation for the specific 
heat at constant volume. 


Temper ature(K) 298.15 400.00 500.00 

C v for n- C 3 F 7 1 [J/moleK] 146 169 186 

C v {or t - C 4 F 9 I [J/moleK] 181 215 240 

Using least squares there results the approximations 

C v = 183.26 exp(0.0014(T - 300)) for t - C 4 F 9 I 

C v = 147.23 exp(0.0012(T - 300)) for n - C 3 F 7 I 

Stimulated Emission Cross-Section 


The stimulated emission cross-section is given by (references 6,7 and 8) 

X 2 Ai 


°ij = 


8tt 




where, gij(v) is the normalized line shape function 

2 


Sfij(i') = 




(31) 


(32) 


The level transitions for the iodine laser are illustrated in the figure 3. The relative 
intensities of these transitions are illustrated in the figure 4. The Einstein coefficients for 
the different lines have the transition rates 


-^34 = 5.0a A 23 —2.3a 

.A33 — 2.1a -^22 =3.0a 

A32 = 0.6a .A21 =2. 4a 


(33) 


in units of sec 1, where a = A/ 7.77 with A = 5.4 ± 2.0sec 1. Using 

1/0 = 1.315246 (10)—* = 2 ' 28 ° 94 (W) U " 2 ' 3 (1 ° )4 GHZ ( 34 > 

the laser emission frequencies from u 0 are given by 


u 34 = i / 0 

v 33 = ^0 + 0.141c 
1/32 = 1/33 + 0.068c 
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i/21 = r/o — 0.427c 
i/22 — ^ 2 1 — 0.026c 
1/23 = ^22 — 0.068c 


(35) 



p l !' 


and consequently the overall stimulated emission cross-section is given by 


a = 


-E 

19 ^ 


^2i 


7 V' 

+ t^E 


^3» 


L 12 H 1 + [2(^-)]] 2 12 1 + [2(^-)] 2 J 

which is based upon the statistical weights of the hyperfine levels (reference 1 ) and 


( 36 ) 


6u = a 0 + a x p (37) 

with <*o = 2.51 (10 ) 8 y/T /300, a x = 1.88 (10 ) 7 \ZtJ3QO where p is the pressure in torr, 
and T is the absolute temperature. Here ao is related to the Doppler line width and ai is 
the pressure broadening coefficient associated with the lasant n — C 3 F 7 J. 

3 n — 

2 Pi / 2 

2 -m 


4 


2 Pz/2 


1 


Figure 3. Level transitions for iodine 


17 



13% 

2-3 


12.3% 

15.2% 

2 - 22-1 


36.5% 16.3% 6.2% 

3-4 3-3 3-2 


^5 


^4 


V3 


WO 


"1 


I/ 2 


1.026 


.068 


.427 


.141 


.068 


Figure 4 Relative intensities for iodine transitions 


The figures 5, 6, 7, 8 contain graphs of the stimulated emission cross-section for various 
pressures. 
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. Figure 5 Stimulated emission cross-section (10)~ 2O cm a 

vs laser emission from Vo> frequency GHz with pressure of 5 torr 

II 


Voigt profile of cross-section in units of 10 
300.0 n Pressure = 5 torr 







The above models were used to compare theoretical output with results from various 
laboratory experiments. It was found that certain of the reaction rate coefficients had 
to be changed in order to make theory agree with experiment. Additional experimental 
results indicated that the purity of the iodide gases was in question and so we felt justified 
in varying some of the reaction rate coefficients. The Appendix A contains a sampling of 
the various theoretical results obtained during the course of this study. 
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Miscellaneous Results 
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n-C3F7I 


t-C4F9I 


PTO = 
OMEG1 = 
COO = 

R1 = 

R2 = 

TM = 
XNRHO = 
CON = 

LC = 

ZOL = 

A = 

R20 = 
FRAC = 
TO = 

RAD = 

VI = 

V2 = 
TTT2 = 
TTT3 = 
TTT4 = 
TTT5 = 
TTT6 = 
CHI1 = 
CHI2 = 
CHI 3 = 
KK1 = 
KK2 = 
KK3 = 
KK4 - 
KK5 * 
KK6 = 
KK7 = 
KK8 = 
KK9 = 
KK10 = 
AAO = 
BBO = 
CC1 = 
CC2 = 
CC3 = 
CC4 = 
CC5 = 
QQ1 = 
QQ2 = 
QQ3 - 
QQ4 - 
QQ5 = 
WTMOL = 


5. 600E+00 
6. 630E+02 
1 • 000E+17 
1 . OOOE+OO 
7 . OOOE-Ol 
3 . 000E-01 
1. 000E+00 
4 . 500E+02 
3 . 300E+01 
7 . 500E+00 
1. OOOE+OO 
0. 000E+00 
3 . 750E-04 
3 . 000E+02 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
1. 000E+18 
1. 000E+18 
1 . 000E+18 
1 . 000E+18 
1 . 000E+18 
1 . 200E-02 
1. 200E-01 
0. OOOE+OO 
1 . 000E-14 
2.300E-11 
2 . OOOE-12 
3 . 000E-16 
1 . 000E-11 
0. OOOE+OO 
3 . 000E-19 
1 . 600E-23 
1. 000E+15 
1. 000E+17 
1 . 472E+02 
1 . 2 00E-03 
1 . 600E-33 
5 . 700E-33 
0. OOOE+OO 
1. OOOE+OO 
8 . 000E-33 
1 . 700E-17 
2 . 890E-11 
3 . 700E-18 
4 . 700E-16 
1. 600E-14 
2 . 960E+02 


PTO - 
OMEG1 = 
COO = 

R1 * 

R2 = 

TM - 
XNRHO - 
CON = 

LC = 

ZOL = 

A = 

R20 = 
FRAC = 
TO - 
RAD = 

VI * 

V2 = 
TTT2 = 
TTT3 = 
TTT4 * 
TTT5 - 
TTT6 = 
CHI1 = 
CHI 2 = 
CHI 3 = 
KK1 = 
KK2 = 
KK3 - 
KK4 = 
KK5 = 
KK6 = 
KK7 - 
KK8 - 
KK9 = 
KK10 = 
AAO * 
BBO = 
CC1 = 
CC2 = 
CC3 = 
CC4 = 
CC5 =* 
QQ1 = 
QQ2 = 
QQ3 - 
QQ4 - 
QQ5 = 
WTMOL* 


9.00E+00 

5.50E+02 

1.40E+19 

1.00E+00 

7.00E-01 

3 . 00E-01 

1.00E+00 

4.50E+02 

3.30E+01 

7 . 50E+00 

1.00E+00 

0.00E+00 

3.75E-04 

3.00E+02 

0.00E+00 

0.00E+00 

0.00E+00 

1.00E+18 

1.00E+18 

1.00E+18 

1.00E+18 

1.00E+18 

1.44E-02 

1.20E-01 

O.OOE+OO 

1.00E-14 

6.00E-12 

3.00E-14 

3.00E-18 

1.00E-11 

0.00E+00 

3 . 00E-19 

1.60E-23 

1.00E+14 

1.00E+16 

1.83E+02 

1.40E-03 

1.60E-33 

5.70E-33 

0.00E+00 

1.00E+00 

8.00E-33 

6.10E-17 

2.89E-11 

3.70E-18 

4.70E-16 

1.60E-14 

3.46E+02 


Table A1 Parameter Values for Figure A1 


n-C3F7I 


p 

Con 

Omega 

Power 

Density 


torr 

S.C. 

cm/sec 

w/cm*cm 

Calculated 

Experimental 

5.6 

450.0 

663.0 

0.550 

0.450 

5.6 

600.0 

663.0 

0.650 

0.573 

6.0 

740.0 

619.0 

0.650 

0.700 

5.8 

925.0 

640.0 

0.770 

0.760 

6.4 

1100.0 

580.0 

0.710 

0.860 




t-C4F9I 



P 

Con 

Omega 

Power 

Density 


9.0 

450.0 

550.0 

2.420 

1.850 

4.5 

600.0 

733.0 

1.610 

2.070 

4.5 

740.0 

733.0 

2.090 

2.290 

4.5- 

925.0 

707.0 

2.700 

2.670 

4.2 

1100.0 

707.0 

2.940 

3.120 


Table B1 Data for Figure A1 


Flow Rate, (M/S) 



Figure A2(a) Parameter fit of t - C 4 F 9 I and n - C 3 F 7 I data (Flow Rate (M/S) Vs Power) 


Theoretical data points 
Experimental data points 
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Figure A2(b) Parameter Fit of t - C+Fvl and n — C3F7J data (Flow Rate (SCCM) Vs Power 


Theoretical data points 
Experimental data points 


n-C3F7I 


t-C4F9I 


PTO = 
OMEG1 - 
COO = 

R1 - . 
R2 = 

TM - 
XNRHO - 
CON = 

LC - 
ZOL - 
A = 

R20 - 
FRAC = 
TO = 

RAD * 

VI = 

V2 = 
TTT2 = 
TTT3 = 
TTT4 = 
TTT5 = 
TTT6 = 
CHI1 = 
CHI 2 - 
CHI 3 ■ 
KK1 = 
KK2 * 
KK3 =* 
KK4 = 
KK5 = 
KK6 = 
KK7 = 
KK8 ■ 
KK9 = 
KK10 = 
AAO - 
BBO = 
CC1 = 
CC2 =» 
CC3 = 
CC4 - 
CC5 = 
QQ1 - 
QQ2 - 
QQ3 = 
QQ4 =* 
QQ5 - 
WTMOL - 


17.000000 
740.100000 

8 . 500000E+17 

1.000000 
6 . 750000E-01 
3 . 250000E-01 

1.000000 

1000.000000 

33.000000 
7.500000 

1.000000 

0 . 000000E+00 
1 . 000000E-03 

300.000000 
0 . 000000E+00 
0 . 000000E+00 
0 . 000000E+00 
1 . 000000E+18 
1 . 000000E+18 
1 . 000000E+18 
1 . 000000E+18 
1 . 000000E+18 
1 . 200000E-02 
1 . 200000E-01 
0 . 000000E+00 
1.000000E-14 
8.800000E-11 
2 . 000000E-12 
3 . 000000E-16 
1 . 000000E-11 
0 . 000000E+00 
3 . 000000E-19 
1 . 600000E-23 
1 . 000000E+15 
1 . 000000E+17 

147.230000 
1 . 200000E-03 
1 . 600000E-33 
5.700000E-33 
0 . 000000E+00 

1.000000 
8 . 000000E-33 
1.700000E-17 
2 . 890000E-11 
3.700000E-18 
4.700000E-16 
1 . 600000E-14 

296.000000 


PTO - 
OMEG1 - 
COO “ 

R1 - 
R2 - 
TM - 
XNRHO * 
CON - 
LC - 
ZOL - 
A - 
R20 - 
FRAC - 
TO * 

RAD - 
VI - 
V2 - 
TTT2 - 
TTT3 - 
TTT4 - 
TTT5 - 
TTT6 * 
CHI1 = 
CHI 2 = 
CHI3 = 
KK1 - 
KK2 - 
KK3 - 
KK4 - 
KK5 ■ 
KK6 ■ 
KK7 =* 
KK8 - 
KK9 ■ 
KK10 - 
AAO - 
BBO * 
CC1 - 
CC2 - 
CC3 - 
CC4 - 
CC5 - 
QQ1 - 
QQ2 - 
QQ3 - 
QQ4 - 
QQ5 - 
WTMOL - 


17.000000 
644.60000 

6 . 150000E+18 

1.000000 
6. 750000E-01 
3 . 250000E-01 

1.00000 

1000.000000 

33.000000 
7.500000 

1.000000 
0.000000E+00 

1. 000000E-03 

300.000000 
0.000000E+00 
0 . 000000E+00 
0 . 000000E+00 
1 . 000000E+18 
1 . 000000E+18 
1 . 000000E+18 
1 . 000000E+18 
1 . 000000E+18 
1.440000E-02 
1.200000E-01 
0 . 000000E+00 
1.000000E-14 
8 . 800000E-11 
3 . 000000E-13 
3 . 000000E-17 
1 . 000000E-11 
0 . 000000E+00 
3 . 000000E-19 
1. 600000E-23 
1.000000E+14 
1. 000000E+16 

183.262400 
1 . 398681E-03 
1.600000E-33 
5.700000E-33 
0 . 000000E+00 

1.000000 
8 . 000000E-33 
6. 100000E-17 
2 . 890000E-11 
3.700000E-18 
4 . 700000E-16 
1.600000E-14 

346.00000 


Table A2 Parameter Values for Figures A2(a)(b) 





n-C3F7I 



p 

Con 

Omega 

Power (W) 

Power (W) 

SCCM 

Torr 

s.c. 

cm/sec 

Calculated 

Experimental 


12.0 

1000 

574.9 

5.28 

4.00 

1700 

14.0 

1000 

695.7 

6.28 

6.00 

2400 

17.0 

1000 

740.1 

♦ 

6.28 

7.20 

3100 


< 

* 

t-C4F9I 



P 

Con 

Omega 

Power (W) 

Power (W) 

SCCM 

Torr 

S.C. 

cm/ sec 

Calculated 

Experimental 


14.0 

1000 

507.3 

10.46 

11.1 

1750 

14.5 

1000 

671.7 

13.85 

12.5 

2400 

17.0 

1000 

644.6 

13.51 

13.8 

2700 




n-C3F7I 



P 

Con 

Omega 

Power (W) 

Power fW) 

SCCM 

Torr 

S.C. 

cm/ sec 

Calculated 

Experimental 


12.0 

1000 

574.9 

5.28 

4.00 

1700 

14.0 

1000 

695.7 

6.28 

6.00 

2400 

17.0 

1000 

740.1 

6.28 

7.20 

3100 




t-C4F9I 



P 

Con 

Omega 

Power (W) 

Power (V) 

SCCM 

Torr 

s.c. 

cm/ sec 

Calculated 

Experimental 


14.0 

1000 

507.3 

10.46 

11.1 

1750 

14.5 

1000 

671.7 

13.85 

12.5 

2400 

17.0 

1000 

644.6 

13.51 

13.8 

2700 


Table B2 Data for Figures A2(a)(b) 
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Figure A3(a) Parameter Fit for » — CzFjI data 


EXPERIMENT 





&PARAM 
PTO ■ 

3.580000 

OMEG1 - 

363.200000 

COO = 

3 . 900000E+18 

R1 - 

1.000000 

R2 - 

8 . 500000E-01 

TM = 

1 . 500000E-01 

XNRHO = 

1.000000 

CON = 

220.000000 

LC = 

15.000000 

ZOL = 

7.500000 

A = 

1.850000 

R20 = 

0 . OOOOOOE+OO 

FRAC * 

4.500000E-04 

TO = 

300.000000 

RAD = 

0 . OOOOOOE+OO 

VI = 

0 . OOOOOOE+OO 

V2 = 

0. OOOOOOE+OO 

TTT2 = 

1 . OOOOOOE+18 

TTT3 = 

1 . 000000E+18 

TTT4 = 

1. 000000E+18 

TTT5 = 

1. OOOOOOE+18 

TTT6 = 

1. OOOOOOE+18 

CHI1 = 

1 . 200000E-02 

CHI2 = 

1 . 200000E-01 

CHI 3 - 

0. OOOOOOE+OO 

KK1 = 

1 . 000000E-14 

KK2 = 

2.300000E-11 

KK3 = 

6 . 500000E-13 

KK4 = 

3 . OOOOOOE-16 

KK5 - 

5 . 000000E-11 

KK6 = 

0. OOOOOOE+OO 

KK7 = 

3 . 000000E-19 

KK8 = 

1.600000E-23 

KK9 = 

1 . 000000E+15 

KK10 = 

1 . 000000E+17 

AAO = 

147.230000 

BBO = 

1.200000E-03 

CC1 = 

1.600000E-33 

CC2 = 

5 . 700000E-33 

CC3 - 

0. OOOOOOE+OO 

CC4 = 

1.000000 

CC5 - 

8.000000E-33 

QQ1 = 

1 . 700000E-17 

QQ2 = 

2 . 890000E-11 

QQ3 = 

3 . 700000E-18 

QQ4 = 

4 . 700000E-16 

QQ5 = 

/ 

CON PC 

1 . 600000E-14 

>WER(W) (THEORY) 

220 

4.42 

450 

9.86 

745 

19.8 

1050 

26.0 

1320 

30.0 


POWER (W) (CALCULATED) 
5.05 
14.7 
23.65 
26.02 
27.2 


Table A3 Data for Figures A3(a)(b) 



Figure A4(a) Simulated t - C+F* Space Laser 


500cm a = 15 cm frae — 0.0005 


Velocity (M/S) 



Velocity (M/S) 
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Figure A4(d) Simulated t — C 4 F 0 Space Laser 


500 cm a = 15 cm frac = 0.0005 





Velocity (M/S) 
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Figure A5(c) Simulated t - C 4 F 9 I Space Laser 
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Figure A5(d) Simulated t - C A F$I Space Laser 


Velocity (M/S) 







Figure A6(b) Simulated t — C+FqI Space Laser 


